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■ Transition to the semiclassical behaviour and the decoherence process for 

Q\ . inhomogeneous perturbations generated from the vacuum state during an in- 

*q | fiationary stage in the early Universe are considered both in the Heisenberg 

and the Schrodinger representations to show explicitly that both approaches 
£h . lead to the same prediction: the equivalence of these quantum perturbations to 

classical perturbations having stochastic Gaussian amplitudes and belonging to 
. ^ ' the quasi- isotropic mode. This equivalence and the decoherence are achieved 

^ ■ once the exponentially small (in terms of the squeezing parameter r&) decay- 

ing mode is neglected. In the quasi-classical limit — ► oo, the perturbation 
mode functions can be made real by a time-independent phase rotation, this is 
shown to be equivalent to a fixed relation between squeezing angle and phase 
for all modes in the squeezed-state formalism. Though the present state of the 
gravitational wave background is not a squeezed quantum state in the rigid 
sense and the squeezing parameters loose their direct meaning due to interac- 
tion with the environment and other processes, the standard predictions for the 
rms values of the perturbations generated during inflation are not affected by 
these mechanisms (at least, for scales of interest in cosmological applications). 
This stochastic background still occupies a small part of phase space. 
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1 Introduction 



A unique and remarkable property of the inflationary scenario of the early Universe 
(irrespective of its concrete realization) is that it opens an exciting possibility to di- 
rectly observe the outcome of a genuine quantum- gravitational effect: generation 
of quasi-classical fluctuations of quantum fields including the gravitational one in 
strong external gravitational fields (in other words, by space-time curvature). Histor- 
ically, this effect was known as "particle creation from the vacuum in a background 
gravitational field", but it is clear now that what can be measured at present are 
not "particles" but rather inhomogeneous fluctuations (perturbations) of the gravita- 
tional field. In non-inflationary cosmological models, the effect of particle creation is 
exceedingly small and does not lead to observable consequences. Just the opposite, 
not only can minimal perturbations of the gravitational potential generated in the 
simplest versions of the inflationary scenario (first quantitatively calculated in [|TJ) 
be sufficient to explain galaxy formation and the large-scale structure in the Uni- 
verse, but also their predicted spectrum (approximately flat, n ~ 1) and statistics 
(Gaussian) have been successfully quantitatively confirmed by the COBE discovery 
of large- angle fluctuations ^f- of the cosmic microwave background temperature 010 
years after the prediction was made. Moreover, in the case of the inflationary scenario 
and in contrast with other cases, the corresponding "pure" quantum-gravitational ef- 
fect - creation of gravitons by background gravitational fields - produces a large relic 
gravitational-wave background with frequencies <C 10 10 Hz in the Universe [||, and it 
is even possible that a modest, but still significant part of the observed large- angle 
fluctuations (at the level proposed in but probably not larger) is due to these 
gravitational waves. Another possibility of observing creation of particles by gravi- 
tational fields might be through the Hawking radiation from primordial black holes 
(PBH) with masses M < 10 15 g, but it follows from direct or indirect observational 
tests (see, e.g. ||) that the number density of such PBH in the Universe is very small 
if they formed at all (and the inflationary scenario typically predicts their complete 
absence) . 

In spite of this definite success of the inflationary scenario and of the quantum 
theory in curved space-time, some confusion still seems to exist in the literature 
regarding how rigid the derivation of the perturbations is (see, e.g. polemics in || 
|7], The key problem here is that though the process of creation from the vacuum 
and the perturbations themselves are purely of a quantum-mechanical nature (at 
least initially), the observed temperature or density fluctuations in the Universe are 
certainly classical. 

Thus a complete derivation should include some mechanism of quantum-to-classical 
transition and decoherence of the perturbations. Connected with this are fundamen- 
tal questions about the wave function of the Universe being pure or mixed and its 
interpretation. An additional complication is the relation between the Heisenberg 
and the Schrodinger representations in quantum mechanics and quantum field the- 
ory. Of course, these representations are completely equivalent and contain the same 
physics, but they use a different language and different parameters for the description 
of a given state of a quantum field. Almost all initial studies of particle creation in 
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cosmology in general 0, W, II], and in the inflationary scenario in particular [|T]. [| 
were performed using the Heisenberg approach. This approach is more convenient for 
the purpose of renormalization ITDl [12[ and a description of the quantum-to-classical 
transition (we shall come back to the latter point below). The use of mode functions 
satisfying a classical wave equation and the Bogolubov transformation for creation 
and annihilation operators is a characteristic feature of this approach. Derivation of 
the perturbations in this approach usually ends up (like in |], f|) by taking these 
mode functions as classical variables with stochastic Gaussian amplitudes but satis- 
fying a certain condition in the regime outside the Hubble radius ("non-decreasing 
modes"). On the other hand what naturally follows from quantum cosmology where 
a homogeneous isotropic background is quantized too, is the Schrodinger represen- 
tation for the wave functions of perturbations (see e.g. |[3|]). This representation is 



usually used also in order to consider the decoherence process. Here one speaks about 
a two-mode squeezed state and describes it with the help of squeezing parameters. 
Though, of course, it is generally well-known that the Bogolubov transformation of 
annihilation and creation operators in the Heisenberg picture just corresponds to the 
evolution of the vacuum state into a squeezed one in the Schrodinger picture (see 



e.g. a mathematical analysis in |14|, [TI| |16|), there still exists a point of discussion 



about how the two approaches are related in the cosmological context and which of 
them is "better". Inspired by the impression that a squeezed state has a non-classical 
behaviour even for large values of the squeezing parameter r, there were even claims 
that the squeezed state formalism gives observable predictions which are superior to 
the usual Heisenberg approach || 0] (actually it does not, see also ||). It is clear 
that a deep understanding of the generation process of perturbations is of utmost 
importance both for further development of the theory of quantum gravity and quan- 
tum cosmology and for observational implications. That is why we reconsider this 
question here. 

We will deal with the simplest case of a quantum real massless scalar field in a 
Friedman-Robertson- Walker (FRW) background because it is sufficiently representa- 
tive for our purposes. First, the equation for the time-dependent part of <fi coincides 
exactly with the equation for the time-dependent part of gravitational waves on a 
FRW background [17 under the condition that the non-diagonal components of the 
matter pressure tensor are zero in the first order of perturbation expansion. The 
latter condition is satisfied e.g. in the case of matter consisting of a mixture of an 
arbitrary number of hydrodynamical components with barotropic equations of state 
Pi — Pi{ti) and scalar fields with arbitrary mutual interactions and minimal coupling 
to the gravitational field. Furthermore, this equation has a structure very similar 
(though generally not completely identical) to the equation satisfied by the gravita- 
tional potential. 

In section 2, we remind the reader of explicit relations between mode functions 
of the field, coefficients of the Bogolubov transformation and squeezing parame- 
ters. Then, in section 3, we consider the quasi-classical limit in the Heisenberg and 
Schrodinger representations in parallel in order to emphasize the fact that a two-mode 
squeezed state for large absolute values of the squeezing parameter |r| is completely 
equivalent to a classical standing wave with a stochastic Gaussian amplitude. The 
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term "standing" means that there is a definite deterministic correlation between k 
and — k modes for each wave vector k. Also, we show how this property is related to 
the fundamental fact that the field modes can be made real by a time-independent 
phase rotation in this limit. All this is illustrated with a specific but very impor- 
tant example, namely that of the de Sitter background. In section 4, we remind the 
physical process in the Universe leading to extreme squeezing (|r| — > oo) and to the 
quantum-to-classical transition - the different behaviour of non-decreasing and de- 
caying modes outside the Hubble radius. Then it follows that quantum-to-classical 
transition and decoherence in the Heisenberg representation (in contrast with the 
Schrodinger one) are achieved simply by omitting an exceedingly small part of the 
field operator (the decaying mode), without any need to consider some interaction 
of the mode with an "environment". This may be called, following J. A. Wheeler's 
favourite way to put it, "decoherence without decoherence". After this omission, it 
becomes unimportant whether the field is in a pure or in a mixed state. As a result, 
we come to the conclusion that the Heisenberg (field mode) approach becomes more 
straightforward in real situations when very small interactions of the field with other 
fields take place. Namely, due to these interactions and the resulting decoherence 
process, the present quantum state of the field is neither a pure squeezed state nor 
even can it be described by a squeezed density matrix, so the use of squeezing param- 
eters looses sense. On the other hand, interactions practically do not change the field 
modes (at least for sufficiently large scales). Thus, all predictions about present-day 
perturbations remain unchanged. A possibility of having nevertheless some "quan- 
tum signature" in the present-day spectrum of perturbations is mentioned. Finally, 
we discuss the point that the omitted, exponentially small part of the field may be 
important for the calculation of the entropy of the perturbations. 

2 Bogolubov transformation and two-mode squeezed 
state 

We give here the essential about quantized fields on a flat FRW background. Let us 
consider a real massless scalar field 0. It is described by the following action S: 



where [i — 0, .., 3, c = h = 1 and the Landau-Lifshitz sign conventions are used. The 
space-time metric has the form 



Let us remind how the dynamics of this system will lead to the appearance of squeezed 
states. We first write down the classical Hamiltonian H in terms of the field y = a<p 
and the conformal time r] = J -Mr. The following result is then obtained 




(1) 



ds 1 = dt 2 - a 2 (t)S ij dx i dx j , 



i, J = 1,2,3. 



(2) 



H 
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= i J d 3 k[p(k)p*(k) + k 2 y(k)y*(k) + a - (y(k)p*(k) + p(k)y*(k))] (3) 

where 

_d£(y,y') , a' 

P = — q ^ — = 1/ 2/ 4 

oy' a 

and a prime stands for derivation with respect to the conformal time. Here the follow- 
ing Fourier transform convention is used: $(k) = (27r) -3 / 2 / $(r)e _jkr cFr for functions 
as well as for operators. In order to avoid too heavy notations, we will often write 
simply y(k), a(k), ... instead of y(k, 77), a(k, rj), ... though the Fourier transforms are 
time-dependent c-functions or time-dependent operators in the Heisenberg represen- 
tation. Due to reality of the field y, we have that y(k) = y*(— k), resp. y^(— k) for 
operators. Therefore, any classical field configuration is completely specified by giving 
the Fourier transforms in half Fourier space. This may be not true in the quantum 
case, and the full Fourier space has to be used if a quantum state of the field is not 
invariant under the reflection k — > — k. However, for the vacuum initial state that 
we will use below, there is no such complication. The Fourier transforms appearing 
in (H) satisfy the equation 



a" 



y"(k) + Ik* - - y(k) = 0. (5) 



When the field y is quantized, the Hamiltonian becomes 
d 3 k t 



r H s k n' 
H = / -— [fc(a(k)a t (k) + a f (k)a(k)) + i-{<J{k)a\-k) - a(k)o(-k))]. (6) 

J Zi (X 

The time-dependent (in the Heisenberg representation) operator a(k) appearing in @ 
is defined as usual: 

o(k) = Uk y(k) + <^P0<)) > (7) 



a(k)+at(-k) JU ^ ; fk ( ^ 



y(k)= v ' \ p(k) = - M /-(a(k)-at(-k) 



so that 



The canonical commutation relations 

[y(x,77), p(x',^)]=^ (3) (x-x') (9) 

imply the following commutation relations 

[y(k,r)) , p\k', V )] = i5^(k-k'), [a(k, V ) , a\k', V )} = 5®(k-k'). (10) 

The last piece in the integrand of (|6]) is responsible for the squeezing. Let us see first 
how it affects the time evolution of the system. We have 



a'(k) \ _ f -ik 2- \ / a(k) 
(aK-k))')~{ a - i j 1 « f (-k) 



(11) 
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Clearly, the general solution of these two coupled equations are 

a(k,r]) 



-k,r]) 



u k (r])a(k, rjo) + v k (r})a^(-k, r? ), 
u* k {vW(-K rjo) +v* k {r])a{k, rj ). 



(12) 



This is just a Bogolubov transformation. Eq. (|T2|) can be interpreted as giving the 
time evolution of the creation and annihilation operators in the Heisenberg represen- 
tation, or as a definition of explicitly time- dependent operators in the Schrodinger 
representation. The commutation relations ( |10|) are preserved under the unitary time 
evolution which yields the constraint 



\Mv)\ 2 - \vk(v)\' 



1. 



(13) 



therefore, allowing the following standard parameterization of the functions u k ,v k : 

e -ie k (v) cos h rfe (r/), 



u k (r]) 



,i(fi k {v)+*Pk(ri)) S i n h r k (7]). 



(14) 



Here r k is the squeezing parameter, tp k is the squeezing angle and 6 k is the phase. 
The relation between these quantities and those introduced in the a — (5 formal- 
10| is the following (if Q k in |L0| is chosen to be equal to k): 



ism 



u k 



a k e- ik \ 



(3ke 



(15) 



Also, the parameters s k} u k} r k introduced in [JTIJ, that proved to be very useful to 
make the adiabatic expansion for large k and to obtain a finite average value of the 
energy-momentum tensor of the quantum field either by the n-wave regularization 



method JT0[ or by the equivalent adiabatic regularization method [I2|, are expressed 
through the squeezing parameters as: 



Uk 
Tk 



\Pk\ 2 = sinh 2 r fc 



a k p* k e~ 2 ^ + al(3 k e 2lkr i = cos2^ fc sinh2r fe , 
i{a k (5* k e- 2ikr > - a* k [3 k e 2ikv ) = - sin 2<p k sinh 2r k 



(16) 



and we use here the notation u k to avoid confusion with the quantity u k used in the 
present paper. Note that s k is equal to the average number of created particles with 
momentum k in the WKB regime (in particular, if a(r]) becomes constant). We don't 
intend here to introduce the notion of particles in a non-WKB regime because it is 
ambiguous and does not lead to interesting results. 

Let us now define the field modes f k {rj) with ?kf k = f kl and Qf k = f k2 , fk(vo) = 
1 / \^2k, we will adopt a similar notation for the quantities p and y, 

y(k) = f k (r])a(k,7] ) + / fc *(r/)a t (-k, r/ ) 



2kf k i(r))y(k,r) ) - \j -f k2 (v)p(K Vo) 



(17) 
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and the momentum modes gk(v)-> 9k{Vo) 
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p(k) = -i[g k (r])a(k,r] ) - g* k (r])a^(-k,r] )] 
f 2 



k 9ki(v)p(^Vo) + v2fc g k2 (r])y(k,7] ). (18) 
The modes fj~ satisfy the Euler-Lagrange equation (|]). Note that 

^ = \[2~k 1 |/ fc | 2 = —(cosh 2r fc + cos 2<^ fe sinh 2r fc ), 

<?* = \fl(u k -v* k )=i(f' k -^f k ). (19) 

Eqs (pT^ , p~9|) give explicitly the relation between the modes which are typically 
used in the Heisenberg approach and the squeezing parameters characteristic for the 
Schrodinger approach. Also, they can be used to obtain the dynamical equations 
satisfied by the squeezing parameters, see Eq. ( |3~9"D below. The Wronskian condition 
for Eq.(|5]), as well as the commutation relations (|H]), yield the following equality 

9k fk + 9lfk = iWkK ~ fkfk) = 1. (20) 

We will be interested, in particular, in the quantum state of the field y defined to be 
vacuum at some time t]q in the following way 

Vk: a(k,»fo)|0,77o) = 0. (21) 

This state corresponds to a Gaussian state and time evolution preserves its Gaussian- 
ity Indeed it follows from (|17])|I1) that in the Heisenberg representation, the time 
independent state |0, t]q)h is an eigenstate of the operator y(k) + ry^T 1 (r])p(k) , namely 

{y(k) + i lk -\ V )p(k)}\0, m ) H = 0, (22) 

where the operators y(k),p(k) as well as the function ^ k depend on time, 

^u* k — v k ^ 1 — « sin 2y9 fc sinh 2r fe 1 .F{k) 

u* k + v k cosh 2r k + cos 2ip k sinh 2r k 2\f k \ 2 % |/ fc | 2 

F(k) = ^su k v k = ^sfk*9k = - sin 2<p k sinh 2r k ■ (23) 

On the other hand, in the Schrodinger representation the time-evolved state |0, 77)5 = 
S\Q,rj ), where 5* is the S'-matrix, satisfies the equation 

Sa(k, Vo )S- 1 \0,r ] ) s = (24) 

or equivalently 

{y(k, m ) + i lk - l (r l )p(k, Vo )}\Q,r ] ) s = . (25) 
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Note the similar structure of Eqs P2]P3| ). In the coordinate Schrodinger represen- 
tation, p(k, T] ) = —id/dy(—)s.,r)o). Hence the state (0,770)5 has a Gaussian wave 
functional in this representation consisting of the product of 

*[y(Mo),y(-Mo)] = iV.exp (- ^^^ {1 - z2F(fc)}) 

= ^«p(-J?^C{l-<2F(fc)}) (26) 

for each pair k, — k where N k is a normalization coefficient. The time dependence 
of \1/ is through f k , F(k), and N k . This structure of the wave functional just re- 
flects the fact that we get a two-mode squeezed state. The corresponding probability 
V[y(k,T]o),y(-k,T]o)] is given by 

P[y(k,77o),y(-k,77 )] ocexp ( J^J^H ■ (27) 

At rj = rjo, we have jkivo) — k or equivalently F(k) = 0, in other words we have a 
minimum uncertainty wave function. 



3 Transition to semiclassical behaviour 

Let us first consider the transition in the Heisenberg approach and take the formal 
limit "H — > 0" keeping the rms amplitude \f k \, when expressed in physical units, 
fixed. Since the right-hand sides of the commutation relation (|9|) and of Eqs ( fl3| , |2~0"|) 
are proportional to h in physical units and do not depend on |/ fc |, they may be 
approximately replaced by in this limit. In other words, \v,k\ ~ \vk\ ^> 1, |/fc| 3> 
l/\^2k, \g k \ 3> \J\ in natural units in the quasi- classical limit. Then it follows from 
Eq. (|20|) with in the right-hand side that /| = c k fk, where c k is a time-independent 
constant. As a result, it is possible to make f k real for all times by a time-independent 
phase rotation, viz. f k — > f k exp(— |argCfc). On the other hand, g k is purely imaginary 
in this limit. 

A further consequence is that all variables y(k) and p(k) become mutually com- 
muting. However, we still cannot ascribe any definite numerical values to them, in 
contrast with coherent states in the quasi-classical limit; there is no Bose condensate. 
The correct way to put it is that field modes become equivalent to stochastic c-number 
functions of time with some probability distribution p(y(k)y(— k)) = p(\y(k)\ 2 ) for 
each pair of modes k, — k in the following sense 

H (0, m \G(y(k))tf(y(k))\0, Vo ) H = J J ^(^^(^^(^(^^^^(k))! 2 (28) 

for any arbitrary function G (y(k)) and k / 0, (for k = the proof is a little bit 
different). Here we assume for simplicity of notation the k spectrum to be discrete, 
this can be achieved e.g. by formally making the T 3 identification of space with three 
topological comoving scales much larger than all scales of interest. By considering 
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average values of arbitrary powers of y(k) and y^(k) and using the Wick theorem, it 
is straightforward to show that the two-dimensional probability distribution p(|y(k)|) 
is Gaussian with the dispersion (yf) + = \fk\ 2 - Indeed, if 

r^O ^IM V \fk\ 2 J 

where the argument k of y is omitted for brevity and fk has been made real, then 

oo oo 

H {0,r] \G(y)G^y)\0,r] ) H = Y.T.^qlfT^ h^M (a(k, Vo ) + a+(-k, m 

m=0 n=0 

oo 

a^k,^) +a(-k,r] )) n \0,r] ) H = ^ m! |g m | 2 /f" v 



'fc 



m=0 



oo oo r r r r 

EE«; d yi dy 2 p(\y\)y m y* n = / / ^^(M)!^)! 2 (30) 

m=0n=0 J J J J 

Thus, in the continuous limit, the equivalent classical stochastic field can be written 
as y(k, rf) = fk(r])e(k) where the quantities e(k) are time-independent 5-correlated 
Gaussian variables with zero average and unit dispersion: (e(k)e*(k')) = ^^•'(k — 
k'), e*(k) = e(— k). The crucial property here is that the time-dependent part fkiji) 
factorizes both in the operator part a(k, r] ) + cfl(— k, r] ) of the field mode y(k, rf) 
and in the stochastic variable e(k) of the equivalent c- number function in the limit 
involved. As a result, after some realization of the stochastic amplitude of the field 
mode has occured, further evolution of the mode is deterministic and is not affected 
by quantum noise. This kind of quantum-to-classical transition is similar to that used 



in the stochastic approach to inflation [18, 19 



We turn now to a description of the same transition in the Schrodinger represen- 
tation. We see from (|26|) that semiclassicality is implied if the following condition is 
satisfied 

|F(*)|>1. (31) 

It is clear from fl3lD that this requires the quantum state to be extremely squeezed, 
namely |r/-| 3> 1. Note that, in this limit, we cannot omit the number 1 appearing 
inside the figure brackets of (^) because it would make the wave function non- 
normalizable. The classicality is to be understood in the following sense: if we assign 
to each point in phase-space (y(k), p(k) ~ SSy(k)) the probability given by p7|), it 
will move with time according to the classical Hamiltonian equations. Writing y(k) = 
|y(k)|e , it follows from (|27|) that \y(k) \ obeys a Rayleigh distribution while the phase 
$ becomes a stochastic variable uniformly distributed (0 < $ < 2tt). Condition fl3T|) 
can also be expressed in the following way 

1 7*2 1 



i > 1 , Ik = 7fci + hk2 ■ (32) 

This can also be seen from the following vacuum expectation values 

(Ay(k,r))Ay\k',r))) = Ay 2 (k)6^\k - k') = |/ fc | 2 ^(k - k') , 
(Ap(k,7/)Ap t (k',»7)) = Ap 2 (k)5^\k-k') = |^| 2 5 (3) (k-k') . (33) 
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where A<3> = $ — ($) and we have further adopted the notation (3>(k, rj)^(k', rj)) = 
^ 2 (k)8^(k — k'), where the quantity $ 2 (/c), the power spectrum of the quantity $ 
depends only on k if the state is invariant under spatial translations and rotations. 
Note further that (y) = (p) = 0. The following identities can then be shown to hold, 
in complete accordance with the wave functional representation : 

Ay 2 {k)Ap 2 {k) = \f k \ 2 \g k \ 2 

= (fki9k2 — fk29ki) 2 + 2 

= i(sin 2 2(f k sinh 2 2r k + 1) = F 2 {k) + - . (34) 

Then the condition of semiclassicality fl3"T| ) corresponds to an uncertainty which is 
much bigger than the minimal one allowed by the rules of quantum mechanics. This 
shows once more that we are dealing in this limit with stochasticity of a classical type. 
It is interesting to note that we also have the following vacuum expectation values 

(A^k^A^k',^)) = JA^(3)(k-k') + 5 (3) (k + k')) , (35) 

(Ay 2 (M)Ay 2 (k',?7)> = JA^(5(3)(k-k')-5 (3) (k + k')) , (36) 

(Ap 1 (k,r ] )A Pl (k',r ] )) = l^!(5(3)( k -k') + 5 (3) (k + k')) , (37) 

(Ap 2 (k,r ] )Ap 2 (k',r ] )) = M(*W(k _ k ') - 5®(k + k')) . (38) 

The equalities (p5| - p8D express the correlation existing between k and — k modes. 
They are in complete agreement with the fact that for any real quantity / with 
f(k) = /i(k) +i/ 2 (k), we have h(k) = A(-k), / 2 (k) = -/ 2 (-k). 

Let us prove now that f k can be made real for \r k \ — > oo using the squeezed 
state formalism. When we write the equation of motions for 9 k and <p k , the following 
equations are obtained 

a' 

r k = — cos 2Lp k , 
a 

i a ' 
if k = —k coth 2r k sin 2ip k , 

a 

@'k = ^ "I tanh r k sin 2(p k . (39) 

CL 

The first two equations are coupled and their solutions can then be substituted in the 
last equation. However, it is interesting to combine the last two equations, yielding 
the following property 

lim (8 k + <p k )' = . (40) 

This means that 9 k + (p k — > 5 k , where 5 k is some constant phase. This in turn implies 
that the field modes f k have the following asymptotic behaviour for r k — » oo 

V2kf k ^e- l5k e rk cosif k , (41) 
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and analogously for the modes g k 




ie~ i5k e rk sincp k . (42) 

The quantity 5 k can be made zero by a time-independent phase rotation in which 
case fk becomes real. It is interesting to note that (|3l"D implies 

lim sin 2<pk sinh 2r k = oo (43) 

even if sin 2(ft — > 0. This is also in accordance with the fact (n^) = sinh 2 rkS^(k — k') 
so that large alone is enough in order to have semi-classicality. 

One more way to study the quantum-to-classical transition is through the Wigner 
function formalism. The Wigner function for the two-mode squeezed state k, — k is 
defined as 

W(y(k),y(-k),p(k),p(-k)) = tA^// dx 1 dx 2 e-^ + ^\y(k)-^\p\y(k)+^) 

(44) 

where p is the density matrix of the state and y\ and y 2 are the real and imaginary 
parts of y(k) in the Schrodinger coordinate represenation, and the same convention 
applies to x(k) and p(k). Using p = with \l/ given by we obtain 



xjky 



W(k, -k) = j^-y 2 J J d Xl dx 2 e -*(Pi*i+P*w)tf* ^(k) - ?QpJ * (y(k) + 

iV^exp (-^) exp (-|/ fe | 2 b(k) - ^(k)l 2 ) -(45) 

W > in this case, so there is no problem with its interpretation as a probability 
distribution in phase space. As explained above, the quantum-to-classical transition 
is achieved by taking the formal limit h — > keeping the physical amplitude \j k \ fixed. 
Then the Wigner function becomes 

w(k, -k) = ivf exp (-^r) * (pi(k) - jj^v^)) 6 (j*0O - {p^ k )) = 



P[y(k), y(-k)} S ( p(k) - T^y(k) .(46) 



This just describes the deterministic motion in phase space of a bunch of trajectories 
with stochastic Gaussian initial amplitude . 

The Wigner function coincides with the probability distribution in phase 
space that was obtained in |2(J for a toy model of an upside-down harmonic oscillator. 
We see that this behaviour of the Wigner function is general in the large-squeezing 



limit and does not rely neither on assumptions made in [2(1, nor on the existence of 



the inflationary stage. In the paper |2T| an objection was raised against this way of 



making the quantum-to-classical transition based on the observation that it is always 
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possible to find new canonical variables y(k),p(k) for which {y }{p ) = \ (even in the 
limit \r k \ — > oo). However, it is clear from the previous discussion that this objection 
is not relevant because the requirement (y 2 )(p 2 ) 3> 1 for all possible canonically 
conjugate variables is not the necessary condition for the semi-classical behaviour. 
The only necessary condition is \r k \ 3> 1, or Eq.flSED below for the physical amplitude 
of gravitational waves. 

In inflationary theories, primordial perturbations are generated by vacuum quan- 
tum fluctuations of a real scalar field where the power spectrum of the quantum 
field fluctuations is given by \fk\ 2 - Let us consider the very important example of a 
massless real scalar field on a (quasi) de Sitter space. In that case we have 

SA = e-*(l-±), ^ = e-*\ 1~<0. (47) 



The modes (|47D give also a very accurate description for slowly varying Hubble pa- 
rameter H, namely when \H\ <C 3H 2 . After some straightforward calculation we 

get 

Uk = e-*^*) cosh r k v k = e i{kv+ ^ sinh r k (48) 

and sinhr/c = — > — oo when krj — > 0. The crucial point in fl4"B|) is that t&n5 k = ^r, 
therefore 5^ will tend to the constant value — |. Hence in the limit krj — > 0, the modes 
f k are purely real up to a constant phase transformation. Let us give for completeness 
the solution for ip k and 9 k 

71 1 1 

cpk = arctan ■ 



4 2 2kr] ' 

= fc^j + arctan — - — . (49) 
2kr] 



We have finally 
for krj — > though sin 2<£>fc — > 0. 



F(ife) = ^ sin 2<^ fc sinh 2r fc ~ (2fc^) _1 -> -oo (50) 



4 Long-wave mode behaviour and decoherence 

The physical mechanism producing the Bogolubov transformation and the extreme 
squeezing in the Universe is, as well known, the expansion of the Universe (not neces- 
sarily inflation) and the existence of the Hubble radius Rh = H~ l = ^j. For modes 
with kR H <C a, i.e. with wavelengths outside the Hubble radius, the general solu- 
tion of Eq. flU) has the following form in terms of the mode functions f k , with y(k) 
expressed through f k using ( |TTD 0: 

f k = d(k)a + C 2 (k)a f V , g k = (^(k^ar?) + i^^- . (51) 
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Ci(k) can be made real by a phase rotation. Then the normalization condition leads 
to the relation 

C^C 2 = - l - . (52) 

The first term in (|5l|) is the quasi-isotropic mode, also called the growing mode. 
It corresponds to a constant value of the field 0. The same behaviour is shared 
by the leading term of scalar (adiabatic) metric perturbations in the synchronous 
gauge for an arbitrary scale factor a{rf) as well as by the gravitational potential $ 
during stages of power-law expansion. The name quasi-isotropic means that this mode 
does not spoil the isotropic expansion of the Universe at early stages, furthermore 
it is contained in the linear expansion of the Lifshits-Khalatnikov quasi-isotropic 
solution P2"| . The second term in (|5*T| ) is the decaying mode. 
In inflationary models, comparing (f47|) with (|51|), we get 



C - Hk C - %k " 2 (53) 

Cl ~7W^ C2 -~V^h (53) 

in agreement with (|52"D where Hk is the value of the slowly varying Hubble parameter 
H at the moment of the first Hubble radius crossing rji (rji < 0, = k~ l ). Here 

we have multiplied fk in (|47|) by —i to make it both real and positive. Then both terms 
in ([51] ) are of the same order at 77 ~ 771. After that, the decaying mode quickly becomes 
exceedingly small. For example, for scales of interest for cosmological applications 
that crossed the Hubble horizon about 60 — 70 e-folds before the end of inflation, the 
rms of the decaying mode is e~ 2rfe ~ 10~ 80 or less from that of the quasi-isotropic 
mode at the end of inflation, and it becomes < 10~ 95 up to the present moment. It is 
clear that we should neglect the decaying mode completely. One more formal reason 
for this is that one cannot keep small terms of the relative order of e~ 2rk as far as 
renormalization contributions to each mode, which are of the order of e~ rk , are not 
taken into account. 

But once the second term in the expression for in (^T|) is omitted, we obtain im- 
mediately decoherence because quantum coherence is described by the correlation (|52"D 
between the non-decaying and decaying modes. Therefore, when working with the 
field modes fk{t]) i n the Heisenberg representation, there is no need to consider any 
interaction with an "environment" and trace over its degrees of freedom in order to 
get decoherence. Moreover, after neglecting the decaying mode, it becomes unimpor- 
tant whether the quantum state of a given mode k (or the Universe as a whole) is 
pure of mixed since the difference between the three following field configurations, 
viz. classical stochastic field with modes /& given by (0), pure squeezed quantum 
state satisfying ( |5Tj , p2| ) and mixed squeezed quantum state with the relation ( |52"D 



understood in the sense of rms values for C*i(k) and C^k), is exponentially small 
(~ e _2rfe ) and disappears after this omission. Summarizing, in the peculiar case of 
quantum cosmological perturbations generated during inflation, the decoherence can 
be obtained without consideration of any concrete decoherence process, that is why 
we may call it "decoherence without decoherence" . 

This property makes the Heisenberg approach more straightforward and conve- 
nient in real physical situations when one takes into account small interactions of 
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the perturbations with matter in the Universe. For example, it is known that the 
interaction of scalar perturbations and gravitational waves with a background mat- 
ter having shear viscosity does not change the quasi-isotropic mode but yields an 
additional exponential decay of the decaying mode. On the other hand, new per- 
turbations belonging to both modes may be generated by local physical processes, 
especially after the second Hubble radius crossing during the radiation- or matter- 
dominated FRW stages. Their amplitude is much more than the fantastically small 
amplitude of the decaying mode remaining from the inflationary stage but much less 
(at least for scales exceeding R eq ~ 15/?r 2 Mpc where h is the present value of the 
Hubble constant in terms of lOOkm/s/Mpc) than the amplitude of the quasi-isotropic 
mode. As a result, the present state of the perturbations is (of course) neither a 
pure squeezed state, nor does it make sense to call the corresponding density matrix 
"squeezed" because one has no reasons to expect to find any direction in phase space 
where the noise (uncertainty) would be less or even comparable to the quantum limit 
| for the vacuum state. A much more adequate description of this state is that it 
consists of the classical stochastic part described by the quasi-isotropic modes /& (the 
first term in Eq.([51]) for k\r]\ <C 1 or its continuation to the regime kr] > 1, see below) 
plus some small noise of indefinite structure. Then the squeezing parameters loose 
their original sense. One may still formally introduce them through the mode func- 
tions using Eqs (|i~4] , |19|) as their definitions (as is done, e.g. in Eq . ([54]) below), but 
they will have little if any relation to the squeezed quantum state in a narrow, rigid 
sense. 

Note that our method of omitting the decaying mode or, equivalently, taking the 
"H — > 0" limit with |/&| being fixed in physical units may be also thought of as coarse 
graining, or smoothing, in phase space. Namely, the radius of the coarse graining 
should be chosen in such a way as to smooth completely the "subfluctuant" variable 
p(k) — r^W(k) (using the terminology of [|l(|) but still without affecting the "super- 
fluctuant" variable y(k) significantly. After that, this radius may be put zero in all 
problems concerning dynamical evolution and stochastic properties of perturbations 
(except for the calculation of their entropy). However, it is clear from the previous 
discussion that final results do not depend on any concrete type of coarse graining as 
far as the leading quasi-isotropic mode remains unaffected by it. 

Let us now consider what happens to the quasi-isotropic mode after the second 
Hubble radius crossing. Then one has fk = Ci(k)f q (k,r]) where f q is the exact 
solution of Eq.(|5p with the asymptotic behaviour f q (k,r]) = a(r]) for krj <C 1, in 
particular, f q (k,r]) = sin kr) during the radiation-dominated regime a = aorj. Note 
that f q = A(k) sin(kr] + for kr] ^> 1 and an arbitrary behaviour of a(rj). Here 
£fc is a constant and A(k) ~ afe) where r] 2 is the moment of second Hubble radius 
crossing, rj 2 ~ k~ x . It is this constant phase or the phase 8^ linearly growing 
with conformal time r\ (see (|55D) that one usually has in mind when saying that the 
gravitational waves generated during inflation have stochastic amplitudes but fixed 
phases. The fact that y(k) becomes zero at some moments of time shows that these 
waves are standing ones. The adiabatic perturbations do not have such an oscillating 
behaviour at the matter-dominated stage. Hence, it is more general to speak not 
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about the fixed phase of perturbations but about the type of the mode (the quasi- 
isotropic one). 

This demonstrates that the loss of quantum coherence does not preclude the exis- 
tence of strong quasi-classical correlations, c.f. comparison of quantum coherence and 
classical correlation in |23| (see also in this connection). If we formally introduce 



the parameters r k} 6 k} <~p k using (|T3]T§), then 



e** = 2kCl[f 2 q + ^((/»') 2 J , n » 1 , 
ip k = -9 k , sm9 k = ) (54) 

for the quasi-isotropic mode. Note the important persisting correlation between <p k 
and 9 k . Deep inside the Hubble radius, 77 ^> r] 2 (k), 

e 2rk ~ 2A;C 2 A 2 = constant , B k = kq + £ k - - . (55) 

It is in this regime that it becomes possible to introduce the number of created parti- 
cles (gravitons) n(k) = n(— k) ~ The condition for the semiclassical behaviour 
of gravitational waves inside the Hubble radius expressed directly in terms of their 
amplitude = —^Jf ~ V~G(f) looks like 

k 3 h 2 g (k) = (A; 3 %(k)/^(k)) ~ e 2r ^ » l f 2 (56) 

where lp is the Planck length and A = 2iTak~ 1 is the wavelength of the perturba- 
tions. The corresponding condition for adiabatic perturbations is the same at the 
radiation-dominated stage inside the Hubble radius, and it is even less restrictive at 
the matter-dominated stage. Thus, if primordial perturbations are measurable at all, 
they are always classical; their quantum origin is reflected in their power spectrum 
and statistics only. It is important to emphasize here that if primordial perturbations 
are quasi-isotropic at present, this does not necessarily imply their quantum origin 
from a squeezed state. It is just the opposite: any classical or quantum process that 
does not spoil the isotropy of the Universe at sufficiently early times results in the 
dominance of the quasi-isotropic mode nowaday. The role of the inflationary scenario 
here is to provide a causal mechanism for the generation of the perturbations, while 
in a non- inflationary cosmology assuming FRW behaviour of the early Universe, 
power spectrum and statistics of the quasi-isotropic mode may be chosen arbitrarily 
by hand. On the contrary, if the present GW background was generated at late times, 
after the second Hubble radius crossing (e.g. by cosmic strings), then both modes 
would be present with equal probability and the phase would be uniformly dis- 
tributed. This, as well known, gives a possibility to discriminate between primordial 
and non-primordial GW. 

Is it possible to verify experimentally the predicted quasi-isotropic character of 
the primordial gravitational wave background on sufficiently large scales, i.e. the 
standing wave behaviour of f k , or that of 9 k in flSSP? A direct LIGO-type experiment 
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is clearly hopeless because it would require a frequency resolution at the level of the 
Hubble constant, i.e. 10~ 18 Hz. It is remarkable, however, that this prediction is 
already proved by observations in the following sense. Let us assume that the present 
gravitational-wave background is of primordial origin and its spectral energy density 
is at least 



10 Vr~J ' Re q <X<Rn 
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~ 10" 14 A 7 < A < R eq (57) 

(e c = 3Hq/8tiG is the critical energy density, R H = Hq 1 = SOOO/i^Mpc, A 7 = 
2nT~ l ~ 1cm), as expected in the simplest versions of the inflationary scenario, or 
larger (the latter is possible for A < 100/i -1 Mpc). This background corresponds to 
the value k 3 ^ 2 h g (k) ~ 10~ 5 at the second Hubble radius crossing. Then it follows from 
the CMB data that this background should be dominated by the quasi-isotropic 
modes oc f q (k) for scales A > 100/i _1 Mpc (corresponding to a multipole number 
I < 40). Furthermore, the success of the primordial nucleosynthesis theory proves 
convincingly that the Universe was isotropic beginning from t ~ Is (T 7 ~ lMeV). 
As a consequence, the present background ( [57] ) should be dominated by the quasi- 
isotropic modes for A > h 2 R eq ~ 15Mpc (£(772) ~ 10 10 s). This region can be further 
expanded using observational limits on the PBH number density [[J. Of course, 
the inflationary scenario predicts the dominance of the quasi-isotropic mode for all 
scales ^> 1cm for the primordial background. As explained above, this quasi-isotropic 
behaviour in itself cannot be interpreted as a proof for the quantum origin of the 
perturbations. However, in more complicated inflationary scenarios it is possible to 
find an observable, in principle, feature of the perturbation spectrum \fk\ 2 (namely the 
absence of zeroes of \fk\ 2 as a function of k) that might be interpreted as a "quantum 



signature" though perturbations themselves are classical nowadays as usual p5 



In addition, the standing wave behavour of fk yields the appearance of a compo- 
nent in the ^ polarization multipoles produced by the primordial GW background 
which is periodic in I. The period of oscillations is Ti = Tt(r]o — i] rec )/i] rec in the 
spatially flat Universe where t]q and 7] rec are the present conformal time and the re- 
combination time respectively. The amplitude of this component is, however, on the 
level ^ ~ 10~ 6 or less, the corresponding effect in the (^") itself is even smaller. 

The background consisting of waves oc f q {k) belonging to the quasi- isotropic mode 
only, despite its stochasticity, occupies a volume of measure zero in phase space. Thus, 
one might expect that its entropy is low, if any. To determine the volume occupied in 
phase space, one has to restore the exponentially small decaying mode, or to consider 
other physical processes leading to the generation of both modes. If the quasi-isotropic 
mode remains unchanged, the resulting effect may be described as the appearance 
of a small stochastic correction to the phases and 9k- There were a number of 
proposals about how to calculate the entropy of the perturbations [^, ^7], |28|], all of 
them leading to the same expression ASk ~ 2r^ (r^ ^> 1) for each mode k appearing 
as a result of coarse graining. It remains unclear, however, if this entropy is the 
minimal possible one that has to be ascribed to the perturbations irrespective of the 
choice of coarse graining. We hope to return to this question elsewhere. 
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